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ABSTRACT 


A  syst «■  is  describe!  which  incorporated  a  Hewlett-Packard 
85  "desk-top"  computer  to  control  a  frequency  synthesizer 
and  read  the  output  of  a  lock-in  analyzer  to  measure, 
display  and  record  the  resonant  frequencies,  amplitudes,  and 
quality  factors  of  several  modes  of  an  acoustical  resonator. 
The  system  is  capable  of  locating,  measuring,  and  tracking 
the  resonant  modes  as  parameters  which  affect  sound  speed 
and  attenuation  are  varied.  An  algorithm  for  rapidly 
fitting  "good  quality"  measured  data  to  a  resonance 
lineshape  is  described  which  determines  quality  factors  to 
precisions  of  better  than  0.1  percent,  amplitudes  to  better 
than  0.01  percent  and  center  frequencies  to  better  than 
O.lppm.  Sample  output  is  provided  for  the  lowest  three 
plane  wave  modes  of  an  air  filled  cylindrical  resonator  in 
the  temperature  range  of  -15  to  25  degrees  Celsius. 
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I.  INTRODUCTION 


A.  BACKGROUND 

Many  applications  in  physios  and  engineering  require  the 
ability  to  determine  tha  canter  frequency  (fo) ,  maximum 
amplitude  (A),  and  qua  liny  factor  (Q)  of  a  resonance.  A  few 
examples  of  this  need  arise  in  the  measurement  of  the  speed 
of  sound,  in  the  raoiprccity  calibration  of  electroacoustic 
transducers  [Ref.  1],  and  in  the  measurement  of  sound 
absorption. 

Prior  to  the  aiirent  of  computers  and  instrumentation 
interfacing,  severe  limitations  were  imposed  upon  the 
experimental  precision  obtainable.  The  variation  of 
amplitude  at  resonance  is  second  order  in  frequency  so  that 
a  small  uncertainty  in  the  maximum  amplitude  results  in  a 
greater  uncertainty  in  the  frequency  corresponding  to  the 
maximum  amplitude.  Also,  the  3  is  calculated  from  the  ratio 
of  a  large  number  (canter  frequency)  to  the  difference  of 
two  large  numbers  which  may  differ  by  a  small  amount  (the  3 
dB  down  frequencies),  resulting  in  an  uncertainty  in  Q  which 
can  be  quite  large. 
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This  project  was  designed  to  produce  an  automated,  high 
precision  system  capable  of  obtaining  and  tracking  the 


center  frequency,  maximum  amplitude  and  2  of  several 
acoustical  resonant  modes  in  any  given  resonator  as  some 
external  parameter  affecting  these  values  is  varied.  The 
approach  was  broken  up  into  four  task  groups. 

1.  Acq  ustica^  Signal  fiegegation  and  Detection 

It  was  necessary  to  determine  which  method  of 
resonance  excitation  and  signal  detection  was  best  suited 
for  automated  data  acquisition  and  analysis.  The  following 
are  the  three  techniques  examined  experimentally: 

a.  swept  frequency  continuous  wave  excitation  with 
trackirq  narrowband  heterodyne  reception  using 
an  HP  3580A  spectrum  analyzer. 

b.  p base- sens  it ive  detection  utilizing  a  two-ohase 
lock-in  analyzer  ana  a  programmable  frequency 
synthesizer  as  a  signal  source. 

c.  broadband  noise  excitation  in  conjunction  with 
an  FFT  analyzer. 


2.  2§i.a  Analysis  Program  D;  vs  loo  men  t 

It  was  necessary  to  determine  a  precise  and 
efficient  computer  algorithm  for  the  fitting  of  the 
resonance  data  to  a  resonance  Liceshape  to  extract  fo.  A, 
and  Q. 
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3.  £lltsi  U£lSllJ&&US& 

This  involved  the  assembly  of  the  required 
components  (computer,  analog-to- digital  converters,  counter, 
etc.)  in  a  custom  rack,  interfacing  and  interconnection  of 
components  to  the  computer,  resonator,  and  each  other,  and 
the  writing  of  the  software  to  allow  the  coaputer/ccntroller 
to  orchestrate  all  of  the  acquisition,  analysis  and  display 
functions. 

4.  System  Test  aqd  SJiluation 

The  final  task  was  to  test  the  system  on  the  lowest 
frequency  plane  wave  modes  of  a  gas  filled,  cylindrical 
resonator  to  determine:  the  precision  with  which  amplitude, 
Q  and  canter  frequency  could  be  measured;  the  tracking 
capabilities;  the  processing  speed;  and  the  utility  of  the 
output. 

B.  CONCLUSIONS  FRDM  TASK  COMPLETION 

Upon  completing  tasks  1  and  2,  it  was  determined  that 
method  b  (frequency  synthesizer/phase-sensitive  detector) 
provided  the  best  noise  rejection  and  resolution,  was 
computer  controllable,  and  resulted  in  a  data  set  which 
simplified  the  impleae ntat ion  of  the  resonance  fit  algorithm 
chosen  from  task  2.  Method  a  of  task  1  was  rejected  due  to 
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the  lack  of  a  suitable 


computer  interface  which  would  allow 
autcaatic  changes  in  frequency  range,  k Iso,  the  sweep  would 
require  catching  the  values  "on  the  fly"  thus  resulting  in 
inaccuracies  which  were  sweep  rate  dependent  and  data  points 
which  were  not  necessarily  equaLly  spaced.  Method  c  was  not 
utilized  due  to  the  inability  to  concentrate  the  energy  in 
the  bands  of  interest,  resulting  in  the  need  for  higher 
power  to  the  resonator  without  an  increase  in  the  signal  to 
noise  ratio  in  the  frequency  domain  of  interest. 

The  completion  of  all  of  these  tasks  resulted  in  the 
system  to  be'  described  in  the  remainder  of  this  thesis.  The 
system  is  capable  of  a  precision  equal  to  or  better  than  one 
part  in  10*  in  frequency,  one  part  in  10,000  in  amplitude, 
and  cne  part  in  1020  for  2. 
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II.  EQUIPMENT  5SI0P  £N£  DESCRIPTION 


A.  EQUIPMENT  SETUP  AND  INTERFAX S 

Figure  1  shows  a  block  diagram  of  the  aquipment  utilized 
in  the  final  system.  The  following  is  a  brief  look  at  the 
overall  function  of  the  system  which  will  be  followed  by  a 
description  of  each  individual  piece  of  equipment. 

The  HP-85  computer  directly  controls  all  of  ohe 
equipment,  with  the  exception  of  the  Lock-in  Analyzer,  via 
the  Hewlett  Packard  Interface  Bus  (HPIB)  .  The  computer 
sends  a  value  for  frequency  aid  amplitude  to  the  HP3325A 
Synthesizer /Function  generator,  which  causes  an  excitation 
of  the  acoustic  resonances  in  tie  resonator.  The  amplitude 
analog  signal  received  by  the  Lock-in  analyzer  is  converted 
to  a  D.C.  voltage  which  is  sent  to  the  HP3497A  Data 
Acquisition/Control  Unit  where  it  is  digitized  prior  to 
being  sent  to  the  computer  via  the  HPIB.  After  processing, 
the  computer  sends  output  to  the  HP2673A  Printer  and  the 
HP7U70A  Graphics  Plotter.  The  3P3456A  reads  the  temperature 
in  degrees  Celsius  from  a  thermistor  mounted  in  the  body  of 
the  resonator  and,  when  interrogated  by  the  HP-85,  sends  the 
current  value  via  the  HPIB.  In  the  test  systea  (Chapter  IV) 
the  HP3325A  outputs  the  required  frequency  and  amplitude 


FIGURE  1  BLOCK  DIAGRAM  OF  SYSTEM 


signal  by  Beans  of  a  power  aaplifier  to  the  elestret 
transducer  at  one  end  of  the  resonator.  An  electret 
transducer  at  the  other  end  of  the  resonator  picks  up  the 
signal  and  sends  it  via  an  amplifier  to  the  E5&S  5204.  The 
EG&G  5204  is  phase* locked  to  tie  HP3325A  and  outputs  a  D.C. 
voltage  which  is  proportional  to  the  Pythagorian  sum  of  the 
amplified  in-phase  and  quadrature  components  of  the  signal. 
This  analog  amplitude  is  sent  to  the  HP3497A  which,  when 
interrogated,  sends  via  the  HPI3,  a  digital  measurement  of 
the  voltage  level.  The  HP-95  also  outputs  values  to  the 
HP2673A  and  HP7470A  to  be  printed  and  plotted  respectively. 

B.  EQ 01  PH 3 NT  DESCRIPTION 

The  Hewlett-Packard  Personal  Computer  (HP-85)  is  an 
8  bit  microcomputer  containing  as  standard  a  16K  byte 
memory,  a  built-in  thermal  printer  and  a  small  CRT.  The 
computer  is  programmed  in  BASIC  computer  language.  As 
utilized  in  this  application,  the  memory  was  expanded  to  32K 
bytes  through  the  U3e  of  an  add-on  16K  byte  memory  module 
(#52903A).  In  addition,  the  HPIB  interface  (I82937A) ,  and 
the  RON  Drawer  (#82  936A)  with  a  plotter/printer  RON 
(#00085-15002)  and  input/outpjt  RON  (#00085-15003)  were 


utilized 


The  HPIB  is  Hswlett-Padcard's  implementation  of  the 
IEEE  Standard  488-197  5  and  is  a  parallel  bos  of  16  active 
signal  lines,  3  data  control  lines  and  5  interface 
management  lines,  it  is  capable  of  interconnecting  up  to  15 
instruments.  The  input/output  and  priater/plotter  ROMs 
further  expanded  the  memory  in  order  to  handle  the 
transformations  necessary  to  transmit  information  via  the 
HPIB  to  the  peripherals. 

2.  filaitaj  XaLtutej  a£345SA 

Alt  tough  the  HP3456A  has  numerous  capabilities,  in 
this  application  it  was  only  utilized  as  a  means  to  obtain  a 
temperature  from  a  thermistor  aid  to  transmit  that  value  via 
the  HPIB  to  the  HP-85  computer.  Its  only  utilization  was  to 
demonstrate  the  system's  ability  to  track  resonances  under 
conditions  of  varying  temperature. 

3.  siaiksaiiic^EiactUs  siittitai  ZllUlk 

The  HP3325A  has  a  frequency  range  from  1  microhertz 
to  20  megahertz  wita  a  resolution  of  up  to  eleven  digits. 
The  peak  to  peak  output  amplitude  can  be  set  from  one 
millivolt  to  ten  volts.  It  is  fully  programmable  via  the 
HPIB  and  was  so  utilized  in  this  application.  The  execution 
time  for  a  frequency  command  via  the  bus  is  7.0  msec  plus 
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2. 8  msec  for  each  digit  plus  2.3  msec  for  a  decieal  plus 


12. 5  esec  for  each  delimiter.  The  execution  tiea  for 
amplitude  is  6.8  esec  plus  2.3  asec  for  each  digit  plus  2.8 
usee  for  a  decieal  plus  13.0  msec  for  a  deliaicer. 

4*  fill!  Sail  SEliiSZA 

The  HP3497A  if  as  utilized  as  a  voltaeter  in  this 
application.  It  received  an  analog  voltage  level  froa  the 
lock-in  analyzer  which  was  proportional  :o  the  full  scale 
aeter  deflection  and  converts!  it  into  HPIB  coapatible 
signals  to  be  transmitted  to  tas  HP-85.  It  was  completely 
program  controlled,  sampling  and  sending  only  when  so 
directed  by  the  HP- 83. 

5.  SG&G  lock-in  3.0 dei  5204 

Because  the  5204  was  not  HPIB  compatible,  desired 
functions  had  to  be  manually  controlled  vice  program 
controlled.  This  included  changes  to  full-scale  sensitivity 
and  tiae  constants.  It  received  a  phase  reference 
synchronization  signal  froa  the  HP3325A  and  the  data  signal 
froa  the  resonator.  Through  vector  aanipulatlon  it  outputs 
the  square  root  of  the  sub  of  the  squares  of  the  inphase  and 
quadrature  coaponants  of  the  aaplified  signal  to  the 
HP3497  A. 
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These  were  ordinary  amplifiers  inserted  in  the 
systea  to  grovida  signal  gain.  The  HP457 A  was  set  for  a 
voltage  gain  of  two  to  drive  tie  resonator  while  the  HP465A 
provided  a  gain  of  +23  dB  for  tie  output  froa  the  resonitor. 

7*  inisllissat.  KaBki.cs  PLkatit  HilkZli 

The  HP2673A  prints  bi-iirect ionally  at  the  rate  of 
120  characters  per  second  with  j  7  Z11  dot  natrix  character 
font  in  a  9  X  15  character  calL  with  twelve  character  sets 
to  choose  froa  and  three  nodes  of  fornattiag.  This  printer 
was  utilized  to  obtain  standard  sized  paper  output  vice  the 
snail  sized  theraal  printer  output  of  the  HP-85.  It  was 
utilized  in  the  fully  program  controlled  node. 

8.  cs  Slot its  I£7470& 

The  HP7473A  was  utilized  to  produce  the  high 
resolution  and  excellent  quality  graphs  Included  in  this 
report  and  suitable  for  publication.  It  is  a  two  pen 
plotter  capable  of  utilizing  up  to  ten  colors  through 
programed  stops.  It  ha3  a  character  plotting  speed  of  up 
to  six  characters  per  second  ia  any  of  five  character  sets 
with  text  written  in  any  direction,  with  or  without 
character  slant,  and  in  varying  sizes.  If  was  utilized  in 


t 

J  t 

the  full;  prog  raised  lode  via  the  HPIB.  One  of  its  nicest 
features  is  the  ability  to  define  the  plotting  area,  thus 
giving  the  operator  the  capability  of  enlarging  or  shrinking 
the  output  size  as  desired,  all  figures  in  this  thesis  were 
produced  on  this  piece  of  eguipiant.  It  i3  to  be  noted  that 
although  these  figures  are  actual  systei  outputs  or  replicas 
of  system  outputs,  they  ware  collected  over  a  period  of  tiae 
which  involved  several  experiments.  Thus,  each  individual 
variable  nay  not  exactly  coincide  in  valua  from  figure  to 
figure  since  paraaaters  such  as  teaperature,  drive 
amplitude,  full  scale  sensitivity,  time  constants,  frequency 
band,  etc.,  all  affect  each  variable  and  were  possibly 
different  for  the  figures  shown. 

9.  Xb2£&isto£ 

The  theraistor  utilized  in  this  set  up  is  an 
HP-0837-0164  resistive  type  with  a  range  of  -80  degrees 
Celsius  to  150  degrees  Celsius  with  an  accuracy  of  plus  or 
ainus  6  degrees  Celsius  for  the  range  of  -75  degrees  Celsius 
to  130  degrees  Celsius.  Because  the  teaperature  was  only 
being  utilized  as  an  indicator  of  change,  it  was  not 
necessary  to  measure  this  quantity  accurately  for  this 
experiment.  Thus,  the  aost  convenient  means  of 
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accomplish!  rg  the  desired  taste  *  as  utilized.  The  thermistor 
has  a  nominal  resistance  of  5333  ohms  at  25  degrees  Celsius 
rising  to  3.684  Hohms  at  *80  degrees  Celsius  and  falling  to 
92.7  ohms  at  150  degrees  Celsius. 

10.  R e sons to£ 

A  detailed  schematic  representation  of  the  resonator 
utilized  to  test  the  system  is  shown  in  Figure  2  [Ref.  2]. 
Although  this  system  oar.  be  utilized  with  any  resonator,  a 
description  of  the  one  actually  used  will  be  given. 

The  resonator  is  a  cylindrical  brass  cavity  with 
diameter  and  length  egual  to  2.54  cm.  The  ends  are  capped 
by  simple  electret  transducers.  A  small  slot  (0.16cm  X 
1.6cm)  at  the  cylinder  miiplane  exists  to  allow  gas  to  enter 
and  leave  the  resonator  (a  helium  recovery  line)  as  utilized 
in  another  application  [Ref.  2]. 

A  permanently  polarized  12  micron  thick  disk  of 
teflon  (the  electret),  aluninizad  on  one  side,  is 
electrostatically  held  against  a  sand  blasted  metal 
electrode,  which  forms  the  electrically  active  element  in 
the  resonator,  while  the  aluminized  side  is  connected  to 


ground 


The  use  of  pra-charged  electret  material  eliminates 


the  need  for  bias  suppliss  ordinarily  present  in 
conventional  capacitive  transducers.  The  cnlored-in  circles 

I 

in  the  figure  are  o-rings  utilized  in  assembling  the 
resonator  tc  provide  a  gas  tight  seal. 
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III.  PROGRAM  iXPLAMA2£ON 


£ 


A  simplified  flow  chart  o f  the  program  is  shown  in 
Figure  3.  A  brief  overview  of  the  flow  will  be  given 
followed  by  a  description  of  the  individual  blocks. 

The  operator  inputs  an  upper  and  lower  frequency,  a 
drive  amplitude  and  a  time  constant.  This  time  constant 
will  be  discussed  later,  but  it  also  needs  to  be  set  or.  the 
lock-in  analyzer  as  well  as  entered  into  the  computer.  This 
then  constitutes  the  input  to  'Search'.  In  'Search',  an 
approximate  center  frequency  and  amplitude  is  displayed  for 
all  resonances  within  the  frequency  band  entered.  The 
desired  modes  are  caosen  by  tne  operator  and  carried  into 
•Sort'.  In  'Sort'  each  of  the  nodes  is  remeasured  and  then 
refined  via  a  five  point  fit  routine  which  results  in  an 
output  of  center  frequency,  maximum  amplitude,  quality 
factor,  and  bandwidth.  A  new  tine  constant  is  determined 
and  all  are  carried  into  'Measure*.  In  'Measure',  a  more 
accurate  center  frequency,  maximum  amplitude  and  Q  are  found 
for  each  mode.  These  then  go  into  'Ravine'  where  they  are 
refined  to  a  high  precision  by  a  modified  Chi-squared 
minimization  algorithm.  The  Last  step  is  'Track'  where  each 
of  the  modes  i3  continually  cheoked  as  well  as  the  resonator 
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temperature  until  a  change  in  any  resonant  amplitude  of  10 
percent  or  a  diffsrer.ee  in  temperature  of  0.2  degrees 
Celsius  is  detected  at  which  time  the  program  loops  back 
into  'Measure'  and  repeats  the  routine. 

It  is  appropriate,  at  this  time,  to  state  that 
throughout  the  program,  any  system  equipment  response  times 
(as  mentioned  in  Chapter  II)  necessary  for  data 
transfer/acquisition  are  allowed  to  elapse  through  the  use 
of  SAIT  statements  prior  to  sampling  any  given  response. 

A.  INPOT 

The  operator  selects  a  start  frequency  and  an  end 
frequency  to  define  the  frequency  spectrum  to  be 
investigated.  Also,  a  drive  amplitude  and  maximum  drive 
amplitude  are  selected,  which  car.  be  any  value  up  to  and 
including  3500  millivolts  (limit  of  the  HP3325A 
Synthesizer /Function  Generator) .  Two  amplitudes  are 
selected  because  later  in  the  program  the  drive  amplitude 
will  automatically  increment  or  decrement  the  value  in  order 
to  keep  the  lock-in  output  in  its  midrange.  The  second 
value  (max.  drive)  is  a  default  value  adopted  should  the 
program  attempt  to  exceed  it  while  incrementing  the  drive 
amplitude.  Finally,  a  time  constant  (tau)  is  selected.  The 
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lock-in  analyzer  divides  the  spectrum  into  256  regions  and 
integrates  over  each.  The  program  calculates  a  maximum  tau 
by 

T  s  256/4*  (fUppef  flowed 

in  order  that  there  is  net  a  space  between  regions  and  that 
the  sensitivity  over  the  region  be  nearly  flat.  This  value 
is  a  maximum  such  that  the  operator  selects  a  value 
available  on  the  lock-in  analyzer  which  is  equal  to  or  less 
than  this  and  enters  the  value  set  on  the  lock-in,  into  the 
computer. 

B.  STARCH 

The  frequency  range  is  divided  into  256  parts  as 
mentioned  above.  Each  of  these  frequencies  is  sent  via  the 
HPIE  from  the  HP-85  to  the  HP3325A  and  a  resonator  amplitude 
value  for  each  of  these  frequencies  is  obtained  by  the  HP-85 
from  the  HP3497A  via  the  HPIB.  There  is  a  programmed  delay 
(WAIT)  of  4  tau  between  the  sending  of  a  frequency  and  the 
samplinq  of  its  amplitude  in  order  to  allow  the  resonator  to 
attain  at  least  98  percent  of  its  steady  state  value  as 
determined  by  the  lock-in  integration  time.  The  choice  of 
256  regions  was  dictated  by  the  horizontal  resolution  of  the 
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HP-85' s  CRT  display.  These  255  values  are  then  plotted  on 
the  CRT  with  the  option  of  a  hardcopy  printout  from  the 
HP7470A  as  in  Figure  4.  This  graphically  displays  where 
each  of  the  resonance  nodes  within  the  specified  frequency 
interval  is  located. 

Fro*  Figure  4,  it  can  be  seen  that  not  only  are  the 
responses  from  three  major  plane  wave  modes  displayed  but 
also  those  of  several  other  (azimuthal,  radial,  mixed) 
modes.  It  is  possible  to  concentrate  on  any  mods  displayed. 
After  the  operator  enters  a  threshold  relative  amplitude 
value,  the  computer  provides  a  printout  on  its  thermal 
printer  listing  all  amplitudes  and  their  associated 
frequencies  which  equal  or  exceed  this  threshold  value  (see 
Table  2  in  Chapter  IV  for  printout).  If  the  operator  is 
satisfied  with  the  output,  then  the  plot  and  the  thermal 
printout  are  utilized  to  input  to  'Sort*  the  approximate 
center  frequency  and  bandwidth  for  each  desired  mode.  If 
the  operator  is  not  satisfied,  the  program  loops  bach  to 
'Input',  new  parameters  are  chosen,  and  the  process  repeats 
itself. 
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C.  SOM 


The  purpose  of  'Sort'  is  to  identify  for  the  computer 
the  nodes  to  be  investigated  and  to  obtain  "reasonable" 
values  of  center  frequency,  maximum  amplitude,  and  quality 
factor  for  each  desired  mode.  "Seasonable",  in  this  case, 
is  defined  as  sufficiently  accurate  to  allow  these  obtained 
values  to  be  automatically  controlled  ir.  *  Measure'.  This  is 
achieved  by  the  operator  inputting  a  center  frequency  and 
bandwidth  for  each  mods  desired  to  be  worked  upon.  As 
configured,  up  to  nine  modes  can  be  selected.  However,  this 
limit  can  be  increased  through  a  program/diaension  change  to 
whatever  number  is  desired.  Mine  was  chosen  as  a  matter  of 
convenience  tc  simplify  memory  management. 

In  this  portion,  the  bandwidth  of  each  mode  is  divided 
into  one  hundred  (an  arbitrary  number  which  gave  sufficient 
resolution)  equal  frequency  steps  which  are  individually 
sent  and  their  respective  amplitudes  collected  in  the  same 
manner  as  in  search. 

Once  again,  a  plot  is  generated  on  the  CRT  with  the 
option  for  a  hard  copy  (Figure  5)  .  If  the  operator  is  not 
satisfied  with  the  oitput,  then  the  program  loops  back  to 
the  start  of  'Sort'  and  new  inputs  are  entered.  If  the 
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PLOT  OF  RELATIVE  AMPLITUDE  VS  FREQUENCY.  MODE 


operator  i3  satisfied,  ths  program  finds  the  maxima 
amplitude  and  its  corresponding  fraquency  (canter 
frequency)  .  This  is  accomplished  by  examining  the  value  of 
the  amplitude  of  each  of  ths  one  hundred  points  and 
selecting  the  one  which  is  greatest.  The  amplitude  value  of 
the  3  dB  down  points  is  then  found  from  this  maximum 
amplitude.  The  program  again  sorts  through  all  of  the  points 
and  finds  the  closest  "bracket in q"  points  to  this  calculated 
3  dB  down  point.  A  linear  interpolation  is  then  utilized  to 
calculate  the  corresponding  frequency  for  the  calculated 
3  dB  down  amplitude  point.  Df  course  this  process  occurs 
twice  since  there  is  both  an  upper  frequency  and  a  lower 
frequency  for  the  3  dB  down  amplitude  value.  This  process 
constitutes  the  "five  pcint  fit". 

From  these  values,  a  quality  factor  is  calculated  which 
is  equal  to  the  center  frequency  divided  by  the  difference 
of  the  upper  and  lower  3  dB  down  frequencies.  This  value, 
as  well  as  the  center  frequency  and  maxiaum  amplitude,  is 
output  to  the  HP-35  thermal  printer  (see  Table  2  of 
Chapter  17  for  printout). 

This  procedure  is  repeated  for  each  of  the  remaining 
modes.  Figures  6  and  7  ar  e  the  plotter  outputs  of  the  other 
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PLOT  OF  RELATIVE  AMPLITUDE  VS  FREQUENCY, 


FIGURE  7 


PLOT  OF  RELATIVE  AMPLITUDE  VS  FREQUENCY, 


two  normal  modes  investigated  in  this  example.  Figure  8  is 
an  example  of  one  possible  reassn  for  desiring  tb  request  a 
new  input  (wrong  renter  frequency).  In  this  case,  it  is 
possible  that  one  of  the  3  dB  down  points  could  not  be  found 
from  the  data  (cut  of  range)  .  From  this  graph,  the  operator 
would  obviously  either  choose  a  new  center  frequency  input 
or  increase  the  bandwidth  such  that  the  3  da  down  point 
would  occur  in  the  data. 

D.  NSW  TAU 

In  this  block,  three  things  occur.  First,  a  new  time 
constant  is  found.  This  time  constant  is  the  acoustic  time 
constant  required  for  each  point  to  attain  63  percent  of  its 
steady  state  value.  This  is  in  analogy  with  the  time 
constant  required  for  a  capacitor  in  an  electrical  circuit 
to  charge  up  to  its  full  value.  The  acoustic  time  constant 
for  each  mode  is  found  by  dividing  the  3  for  the  mode  b v  PI 
times  the  center  frequency  of  the  mode.  Because  of  the 
manner  in  which  Q  is  found,  this  equates  to  dividing  one  by 
PI  times  the  difference  of  the  3  dB  down  upper  and  lower 
frequencies.  After  calculating  the  time  constant  required 
for  each  mode,  the  new  system  time  constant  selected  is  the 
greatest  one  required,  thus  insuring  that  the  worst  case  is 
adequately  covered. 
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FIGURE  8 


IT  OF  RELATIVE  AMPLITUDE  VS  FREQUENCY. 


Second,  the  background  noise  level  is  determined  by  the 
following  method.  The  synthesizer  outputs  a  frequency  off 
resonance  into  the  resonator  and  the  lock-in  analyzer  is 
sampled  one  hundred  times.  The  mean  and  the  standard 
deviation  of  this  set  of  measurements  is  calculated.  The 
standard  deviation  is  equated  to  the  noise  of  the  system. 
It  is  utilized  later  in  determining  the  signal  to  noise 
(SNR)  ratio  for  each  mode.  It  is  also  utilized  in  the 
calculation  of  the  third  item  of  this  blook  which  will  be 
described  next. 

The  final  (third)  item  is  the  examination  of  each  mode 
to  determine  the  maximum  number  of  points  which  can  be  taken 
without  having  the  possibility  of  the  amplitude  value  of  two 
adjacent  points  being  equal  or  in  the  worst  case  inverting 
position  due  to  associated  uncertainties  or  random  system 
noise.  If  the  frequency  step  size  (  A  f)  is  large  enough, 
then 

dA  *  Amo r  -  k  (  f  o+Af ) 

can  be  chosen  so  that  the  inversion  will  not  occur. 
Assuming  a  Rayleigh  distribution  lineshape  then 
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A(f)  =  A/Q 


Replacing  f  with  (fc  *  Af)  and  e  =  Af/zo,  the  result 
equation  is 


r 


AA  *  A  / 


1 _ 

l(-«Q>2+  Z],/j 


For  the  system  descriaed. 


24  paints  are  normally  taker,  i 


frequency  range  of  2  BW  -here  2  =  fo/3<i.  Therefore, 


following  eguation  results: 


2€Q  =  2  Af/3W  =  1/6  <  1 


and,  that  quantity  squared  is  C<1.  Expanding  and  divii 
fay  A 
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If  A  is  chosen  to  be  some  number  of  standard  deviations,  <*■  , 
based  upon  the  noise  measurement,  then 


Ng  N 

a  *  sm 


where  9/\  is  1/S33  (s ignal-to-noise  ratio).  Mow,  Af  is 
found  to  be 


A  f 


>J~ _  _h_ 

»  2  (SNR)  Q 


If  I  is  the  number  cf  points  to  os  taitsn  in  2  3W  then 


8(  SNR) 

N 

If  the  probability  of  an  inversion,  as  described  above,  is 
set  at  0.1  percent,  then  N  =  3.3  based  upon  a  standard 
distribution  curve.  It  is  worth  noting  at  this  point,  that 
the  probability  cf  inversion  away  from  maximum  is  everywhere 
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lower  until  the  signal  level  is  or.  the  order  of  the 
background  noise.  Utilizing  the  0.1  percent  criteria, 
I  <  *  1. 55  *  SQB  (SNS|  .  This  number  (I)  is  output  to  the  HP-85 
thermal  printer.  The  program  then  defaults  to  twenty- four 
points  as  a  minimum.  This  value  was  based  upon  a 
calculation  from  the  Nyquist  sampling  rate.  The  minimum 
from  the  Nyquist  sampling  rate  (in  integer  value)  would  be 
three;  but  it  is  generally  agreed  upon  in  engineering 
applications  that  a  factor  of  three  or  four  is  required  in 
order  to  obtain  a  satisfactory  reproduction  of  the  signal. 
At  this  point  the  bandwidth  was  doubled  to  insure  inclusion 
of  the  3  aB  down  points  ar.d  thus  the  minimum  number  of 
points  obtained  equalled  24.  For  this  choice  the 

requirement  that  an  inversion  ins  to  noise  has  a  probability 
of  0.1  percent  or  less  requires  the  SNR  to  equal  or  exceed 
237  (♦47.5  dB)  . 

The  system  has  worked  wnen  the  indicated  number  cf 
points  was  as  low  as  10  (still  defaulting  to  a  choice  of  24) 
but  occasionally  gave  faulty  values  due  to  the  reasons 
discussed  above.  Although  it  lid  not  always  happen,  it  is 
recommended  that  whenever  the  indicated  number  of  poiirts 
fails  below  24,  either  the  SNR  is  improved  (increase  drive 
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amplitude,  etc.)  or  the  mode  is  net  suitable  for  further 
investigation.  A  sample  of  tie  thermal  output  format  for 
tau,  noise  and  numoar  of  points,  is  shown  in  Table  2, 
Chapter  IV. 

E.  MEASURE 

Each  resonance  mode  has  a  frequency  band  calculated  for 
it  ir.  'Measure1.  This  band  is  twice  the  band  defined  bv  the 

j  * 

3  a3  down  bandwidti.  Also,  the  relative  value  of  the 
amplitude  for  the  center  frequency  is  evaluated  to  ascertain 
the  need  for  adjustment  in  the  drive  amplitude.  This  is 
achieved  by  driving  the  resonator  at  the  canter  frequency  of 
each  mode  and  sampling  the  amplitude  value.  If  it  is 
between  0.3  and  0.95  as  measured  on  the  Data  Acquisiton 
System,  then  no  adjustment  is  necessary.  If  it  is  not 
within  these  limits,  it  is  either  doubled  cr  halved 
depending  upon  whether  it  was  above  or  below  the  limits. 
This  new  value  is  again  tested  in  the  above  manner  until  it 
is  either  within  the  liiits  oc  has  exceeded  the  maximum 
drive  amplitude  which  was  sat  back  in  'Input'.  If  this 

happens  (exceeding)  ,  than  the  system  defaults  to  the  set 
value  of  maximum  drive  amplitude.  Thi3  adjustment  is 
necessary  because  the  full  scale  sensitivity  setting  on  the 


lock-in  analyzer  is  a  manual  rather  than  computer  controlled 
function. 

In  'Measure'  the  resonator  amplitudes  are  obtained  in 
the  same  manner  (frequency  sent--amplitude  measured)  as  in 
•Sort',  only  this  time,  twenty-four  frequencies  vice  one 
hundred  are  utilized  and  the  'WAIT'  becomes  12  tau  which 
allows  the  system  to  reach  within  6  ppm  of  the  steady  state 
value  to  fully  utilize  the  5  1/2  digit  A-to-u  conversion  of 

the  H?  4397A  (as  explained  previously  at  the  end  of  'NEW 
TAJ').  Once  the  data  are  collected,  the  following  seven 
pcir.t  fit  is  applied.  The  maximum  amplitude  point  and  the 
closest  adjacent  points  on  either  side  of  it  are  used  to 
define  a  parabolic  curve.  The  maximum  amplitude  of  this 
parabola  and  the  corresponding  center  frequency  are  found  by 
setting  the  first  derivative  of  the  amplitude  with  respect 
to  frequency  to  zero.  Solution  of  the  resulting  equations 
and  subsequent  insertion  of  the  solution  into  the  original 
equation  yields  a  unique  maximum  amplitude  and  center 
frequency.  This  amplitude  is  then  utilized  in  the 


calculation  of  the 

3  dB 
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amplitudes. 
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corresponding  to  the 
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then  found  by 

linear  interpolation 

in 
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same  manner 

as  previously 
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described  in  'Sort'.  A  new  quality  factor  is  calculated 
utilizing  these  new  frequencies  and  all  three  of  these 
values;  2,  maximum  amplitude  and  center  frequency  are  sent 
into  ' Ravine'. 

P.  RAVI  HE 

The  values  for  maximum  amplitude,  center  frequency  and 
quality  factor  determined  by  the  seven  point  fit  described 
in  the  previous  section  ('Measure')  are  substantially  more 
precise  than  those  possible  using  conventional  techniques 
employing  a  wave  analyzer  or  spectrum  analyzer  under  manual 
control.  However,  the  high  speed  numerical  processing 

ability  of  the  digital  computer  allows  a  greatly  expanded 
treatment  of  the  data  for  only  a  minimal  investment  in 
actual  processing  time. 

Implicit  in  the  assignment  of  A  (amplitude)  ,  fo 
(frequency)  ,  and  Q  (quality  factor)  as  the  three  parameters 
which  characterize  a  resonance  line  shape  (as  is  done 
throughout  this  thesis)  is  the  assumption  that  the 
amplitude,  as  a  function  of  frequency,  near  resonance  can  be 
described  by: 
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&  ( fn)  *  A/Q  [(  fn/f  o-f  o/f  q)  2  ♦  (1/Q)21* 

where  :  A(fn)  is  the  resonance  lineshape 

A  is  the  calculated  maximum  aaplituda 
Q  is  the  calculated  quality  factor 
fo  is  the  calculated  center  frequency 
fn  is  the  measured  frequency  (24  per  mode) 

Using  all  twenty-four  points,  it  is  possible  to  obtain  the 
sum  of  the  squares  of  deviations  of  the  data  from  the 
resonance  lineshape  defined  by  a  particular  choice  of  A,  fo 
and  Q  as: 

24 

S  2  (A ,  f  o  ,Q)  «  ^[Ar.  -  A  ( f  n ) ]  2 
n*l 

where  An  is  the  measured  amplitude  of  the  n(th)  data  point 
and  A(fn)  is  as  above. 

Conventional  data  analysis  then  accepts  the  "best 
choice"  for  A,  Q,  and  fo  as  that  which  causes  S2  to  be  a 
minimum  [Bef.  3].  This  minimization  for  a  linear  function 
is  a  straight  forward  algebraic  process.  However,  an 
analytic  least  square  fit  to  a  non-linear  function  requires 
matrix  inversions.  This,  in  turn,  not  only  requires  large 
blocks  of  memory  and  computer  time,  but  also  can  easily  lead 
to  errors  due  tc  internal  round-off  or  truncation  of  digits 
within  the  computer  [Ref.  4].  For  these  reasons,  the 
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following  method  was  devised  and  utilized  for  each  of  the 
resonant  modes  being  investigated.  Briefly,  S2  was  computed 
for  small  variations  in  A,  fo  and  Q  as  well  as  for  the 
originally  computed  values  for  these  terms  from  the  seven 
point  fit.  These  S2  terns  were  then  utilized  in  a  parabolic 
fit  in  an  attempt  to  find  those  values  of  A,  fo,  and  Q  which 
wculd  cause  S2  to  be  a  minimum.  Because  the  results  of  the 
seven  point  fit  were  already  quite  good,  the  procedure, 
briefly  described  above  and  expanded  upon  in  the  remainder 
of  this  section,  converges  very  rapidly  to  an  astonishingly 
high  degree  of  precision.  This  precision  will  be  analyzed 
and  discussed  in  detail  in  Chapter  V. 

For  the  remainder  of  this  paper,  the  sum  of  the  squares 
of  deviations  (S2)  will  be  referred  to  as  the  C-G  value, 
which  was  chosen  for  historical  reasons. 

One  at  a  time,  each  of  the  three  parameters  (A,  fo,  Q) 
in  the  resonance  Lineshape  (A(fn))  are  varied.  The 

following  is  an  example  to  clarify  this  process. 

First  the  C-G  term  for  the  calculated  3  is  found.  Now  Q 
is  varied  by  +0.5  percent  and  a  second  c-3  term  is  found. 
Finally,  3  is  varied  by  -3.5  percent  and  a  third  C-G  term  is 
found.  These  three  C-G  terms  are  used  to  define  a  parabolic 
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curve  (Pigure  9)  and  a  new  value  for  Q  which  minimizes  the 
C-G  tera  is  found. 

This  new  value  of  2  replaces  that  found  in  the  seven 
point  fit  and  becoaes  the  new  value  to  be  inserted  in  the 
resonance  lineshape.  This  tins,  however,  A  is  the  variable 
to  be  adjusted.  Just  as  occurred  with  2,  C-G  terms  are 
found  for  A,  A  ♦  3.2.IA,  and  A  -  0.25SA.  Again,  these  C-S 
terms  are  utilized  in  a  paraoolic  curve  minimization  and  a 
new  value  for  A  is  obtained  which  now  replaces  the  A  found 
in  the  seven  point  fit.  Finally  *:hese  new  values  of  2  and  A 
are  utilized  to  determine  a  new  fo  in  the  same  manner.  This 
completes  one  iteration  of  the  process  and  thus  provides 
these  new  values  of  2,  fo,  aid  A  for  the  next  iteration. 


Each 

iterat 

ion  refines  the  preci 

sior.  of  the  values 

of  2, 

fo. 

and 

A  and 

provides  this  new 

set  of  values  for 

t  he 

next 

iteration. 

This  process  is 

repeated  until 

the 

next 

iteration  contributes  a  negligibly  small  variation  in  the 

C-G  value  CBef.  3]. 

In  the  actual  program  utilized,  the  terms  were  varied  in 
exactly  the  order  described  because  in  the  relative 
precision  of  the  tens,  2  was  tie  least  precisely  known  term 
(historically,  this  is  always  true).  Since  this  was  the 
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FIGURE  9 


case,  varying  it  first  resulted  in  the  greatest  effect 
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next  lowest  precision  term  was  A,  so  it  was  obviously  the 
next  choice  to  be  varied.  Finally,  fo,  which  was  actually 
well  known,  was  varied.  Beoaise  of  the  relative  precision 
in  Q,  A,  and  fo  the  variations  were  arbitrarily  chosen  at 
the  indicated  values  (0.5X  of  Q,  0.2"  of  A,  and  0.305  5  of 
3H) ,  as  a  natter  of  experience  since  these  values  reflected 
a  sufficient  variation  for  the  'RAVINS'  technique  to  be 
effective.  The  better  the  approximations  for  Q,  A,  ani  fo 
put  into  'RAVINE*,  the  smaller  these  values  can  be  made. 

Upon  completion  cf  the  final  iteration  (program 
dependent)  the  three  newest  values  for  3,  fo,  and  A  are 
utilized  to  find  a  final  C-G  tern  which  is  then  printed  on 
the  system  printer  aid  stored  on  a  data  tape.  This  c-t  term 
is  related  to  Chi-squared  ani  may  therefore  be  used  to 
determine  the  "goodness  of  fit"  [Hef.  3].  At  this  time, 
other  parameters  and  variables  are  also  printed  on  the 
system  printer  and  stored  on  the  data  tape  (time, 
temperature,  pressurs,  mole  nuaoer,  fo.  A,  Q,  SNR,  drive 
amplitude,  sequential  data  tape  item  number,  and  C-G  Ravine 
value)  .  These  will  be  discussed  at  the  end  of  Chapter  IV 
and  a  sample  output  will  be  shown  (Table  1)  . 
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If  there  are  additional  moles  to  be  evaluated,  then  a 
loop  is  made  back  to  'Measure'  and  the  next  mode  is 
evaluated  in  the  saae  manner  as  described  above.  Open 
completion  of  the  final  mode  to  be  evaluated,  the  program 
shifts  into  'Track'. 


G.  TRACK 


In 

this  portion  the 

currently 

mea 

surer 

amplit 

uda 

of 

each 

mode  is 

compared  with 

the  measu: 

red 

value 

cbtai 

ned 

iu 

ring 

'Heasur 

e' . 

If  the 

value  dir: 

:  ers 

by 

♦  /- 

10 

per 

cent 

(indicating  a 

change  in 

some  axt= 

?  r  r.  a 

1  para 

meter) 

t  * 

:hen 

the 

program 

loops 

back  to  ' 

Measure'  to  fi 

nd  the 

new  Q 

#  A  i 

,  an 

d  fo 

for  each  mode.  Also,  the  temperature  of  the  resonator  is 
measured  and  if  a  difference  of  0.20  degrees  Celsius  is 
found  then  the  system  returns  to  'Measure'.  If  neither  the 
amplitude  nor  the  temperature  has  met  the  criteria  necessary 
to  return  the  system  to  'Measure',  then  the  Track'  section 
continually  loops  tarough  itself  until  either  cf  these 
criteria  is  met.  The  system  remains  in  this  endless  loop  of 
'Measure',  'Ravine',  and  'Track'  until  the  operator  causes 
the  system  to  halt,  the  data  overflows  the  available 
storage  space  on  the  data  tape  (850  aode  sr.tri=s)  or  a 
catastrophic  system  failure  occurs  (loss  of  power,  etc.). 
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The  preferable  method  cf  stoppage  is  to  utilize  a  'PAUSE* 
command  from  the  computer  keyboard  sites  this  enaoies 
continuation  or  program  restart  as  desirad.  It  is  to  be 

noted  that  any  keyboard  action  will  result  in  a  halt. 

This,  as  well  as  any  of  the  other  mentioned  halts  (except 
'PAUSE')  may/will  rejaire  re-loaiir.g  the  program  into  memory 
prior  to  re-running  the  system. 
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IV.  SAMPLE  EXPERIMENT 


In  order 

to 

test 

the  sysc 

an 

'  s  abili 

ty  to 

measure 
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ed  quantit 
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e  and 
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cess 

was 

halted . 

The 

pci 
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e  data 

tape 

was 
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o  pr 
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e  Figures 

10 

through  14  c:i 
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stain 

plot 

ter  to  demonstrate 

the  utili 

ty 

of  the 

system. 

Figure  10 

show 

s  a 

point  by 

po 

int  plot 

of  tea 

per atur 

a  in 

degrees  Kelvin 

versus 

time  in  s 

aconds. 

After  t 

iae  10, 

000, 

the  curve  begins  to  appear  as  a  solid  line  because  the 
individual  points  become  too  closely  grouped  together  to  bs 
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visually  distinguished  in  the  plat.  This  likewise  occurs  at 
varying  places  in  the  other  figures.  also,  in  Figure  10, 
there  is  a  slight  dip  at  the  beginning  of  the  curve  because 
the  system  was  activated  as  the  temperature  was  still 
slightly  decreasing  vice  stable  as  descriosd  aocve. 

Figure  11  shows  a  plot  of  the  log  of  the  difference 
between  a  reference  temperature  (at  infinite  time)  and  each 
individually  measured  temperature,  versus  time.  This  plots 
as  a  nearly  straight  line  whose  slope  is  a  measure  of  the 
exponential  thermal  eguilibrium  time. 

Figures  12,  13,  and  14  show  plots  of  the  square  of  the 
normalized  frequencies  of  eaci  mode  (fo/M*)  2  versus  the 
absolute  temperature.  For  an  ideal  gas  these  plots  should 
be  straight  lines.  In  a  perfect  resonator,  each  successive 
resonant  mode  would  oe  an  exact  multiple  of  the  first  and 
the  normalized  frequencies  would  lie  identically  upon  one 
another.  In  this  experiment  mode  two  is  parallel  to  mode 
one  but  slightly  higaer  in  normalized  value  while  mode  three 
lies  diagonally  in-between  mode  one  and  mode  two.  Because 
the  exact  physical  properties  nf  .he  resonator  itself  were 
not  being  investigated,  this  was  not  persued  further.  The 
slight  curvature  of  the  lines  is  probably  due  to  the  fact 
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FIGURE 
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that  the  thermistor  does  not  indicate  absolute  temperature 
exactly.  The  slop?  of  the  lines  and  the  length  of  the 
resonator  could  be  utilized  to  determine  the  constants  in 
the  ideal  gas  equation  of  state. 

These  figures  (13  through  1*1  demonstrate  a  few  of  the 
possible  means  of  displaying  the  data  stored  on  the  tape. 
It  is  intended  that  future  users  will  determine  what  plots 
are  needed  and  generate  them  acoordir.gly.  The  program  used 
to  produce  cue  of  these  plots  from  the  data  tape  is  included 
in  Appendix  3. 

To  enumerate  what  data  is  availabJ-  on  the  tape,  a 
partial  reproduction  (dot  matrix  output  is  not  permitted  in 
this  thesis)  cf  th?  system  printer  output  for  another 
experiment  is  shown  in  Table  1.  The  time  shown  is  the 
average  time  at  which  the  measurements  were  taken.  The 
temperature  is  the  average  at  which  time  the  measurements 
were  taken  (in  degr?es  Celsiusi  .  The  pressure  is  set  to 
zero  because  a  varyiag  pressure  was  not  to  be  evaluated,  and 
thus,  a  means  of  detecting  pressure  was  not  included  in  the 
equipment  setup.  The  next  entry  in  the  line  is  the  mode 
number  for  which  the  data  applies,  followed  by  the  final 
ravined  values  for  fo.  A,  and  3.  The  SNH  is  obtained  from 
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the  division  of  A  by  the  value  found  for  the  noise  (as 
explained  in  Chapter  III) •  Fhe  next  value  on  a  normal 
printout  would  be  the  drive  amplitude  of  the  system  followed 
by  the  sequentially  increasing  amber  of  the  data  tape  item. 
These  two  items  were  deleted  from  this  figure  in  order  to 
conform  to  thesis  requirements  without  requiring  photo 
reduction.  The  whole  line  constitutes  one  data  item  so  far 
as  the  tape  limit  of  350  data  items  are  concerned.  However, 
each  item  in  the  line  is  an  actual  separate  entry  within 
each  of  the  data  iteas.  The  final  entry  is  the  C-G  'Ravine' 
value. 

Table  2  is  a  reproduction  ri  the  cutout  obtained  from 
the  HP-35  thermal  printer  during  the  course  of  an 
experiment.  For  farther  amplification  of  each  of  the 
entries,  the  reader  is  directed  to  Chapter  III,  sections  B, 
C,  and  D.  The  system  printer  was  utilized  for  this  figure 
again  because  dot  matrix  outputs  are  not  permitted  in  this 
thesis. 
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TABLE  2  SAMPLE  OUTPUT  OF  HP-85  THERMAL  PRINTER 


AMPLITUDE  IN  VOLTS 

FREQUENCY 

.01082075 

6817.00 

.00829425 

6896.00 

.01050175 

13789.00 

.00820725 

13848.00 

.00828725 

20405.00 

.01466500 

20484.00 

.02054900 

20563.00 

.01539075 

20642.00 

.00809400 

20721.00 

MODE  l  CENTER  FREQ  IS  0830  AND  AMP  IS  .0181355  Q  IS  190.123117729 
MODE  2  CENTER  FREQ  IS  13750  AND  AMP  IS  .0127195  Q  IS  135.493980730 
MODE  3  CENTER  FREQ  IS  20540  AND  AMP  IS  .  02531475  Q  IS  194. 224445807 

TIME  CONST  >  »  & 8B0824885U  NEW  TIME  CONST  IS  10 
Tha  man  Is  .000378561  Tha  standard  dav  la  9l  7975285484E-7 


NO.  OF  POINTS  FOR  MODE  1  IS  212 
NO.  OF  POINTS  FOR  MODE  2  IS  177 
NO.  OF  POINTS  FOR  MODE  3  IS  250 


NO.  CHOSEN  IS  24 
NO.  CHOSEN  IS  24 
NO.  CHOSEN  IS  24 


61 


7.  muimaa  Aaa  sisaaflEaBitigiis 


This  section  describes  the  method  utilized  to  deter mine 
the  "goodness  of  fit"  and  the  resulting  precision.  The 
pregrain  previously  described  was  utilized  in  a  steady  state 
condition  (no  variation  of  external  parame+ers) .  Also, 
'Ravine'  was  programs  ed  for  fifteen  iterations  per  node, 
with  a  printout  of  the  3,  A,  and  fo  values  for  each 
iteration.  In  this  lamer  the  change  per  iteration  could  be 
evaluated  for  each  of  the  variables.  As  stated  previously, 
'Ravine'  is  completed  (C-G  nir.imrzed)  when  successive 
iterations  produce  a  negligible  variation  in  the  value 
obtained  [  Sef .  3  ]«  ?  or  a  given  sample  run,  the  variation 
between  the  fourteenth  and  fifteenth  iteration  ((15-14) /15) 
was  approximately  one  in  one  miLlior.  (mode  1)  to  one  in  ten 
million  (mode  2)  for  3,  one  in  two  million  (mode  1)  to  one 
in  fifteen  million  (mode  3)  for  A,  and  one  in  six  billion 
(mode  1)  to  one  in  one  hundred  thiry-sevea  billion  (mode  2). 
It  is  evident  that  the  change  per  iteration  is  definitely 
small  and  thus  the  fifteenth  iteration  values  were  utilized 
as  a  base  line  for  tna  calculation  cf  precision  to  follow. 
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TABLE  3 


COMPARISON  OF  PRECISION  OF  ITERATIONS 


3rd/ 15th 
5th/15th 


3rd/ 15th 
5th/15th 


3rd/15th 
5th/ l 5th 


1 


! 


MODE  l 
Q 

1. 04E-3 
3.85E-4 

A 

4.8BE-4 
1. 72E-4 

f 

l.  20E-7 
4.45E-8 


MODE  2 
Q 

8.94E-5 
3. 31E-5 

A 

2.85E-5 
1. 05E-5 

F 

3.62E-9 
l.  34E-9 


MODE  3 
Q 

1. 25E-4 
4.64E-5 

A 

3.  996-5 
l.  47E-5 

F 

1.54E-9 
5.  71E-9 


2ach  C-G  ravin®  value  for  the  third  and  fifth  iteration 
was  compared  to  the  base  line  to  determine  the  precision  of 
the  system.  This  as  accomplished  by  talcing  the  ratio  of 
the  difference  of  the  base  line  and  the  iteration  to  the 
base  line  value.  Table  3  she**  the  approximate  (rounded) 
values  obtained  for  each  by  mode.  T.nis  represents  a 
best/worst  precision  of  one  in  ten  thousar.d/one  in  one 
thousand  for  Q,  of  one  in  thirty  thousand/one  in  two 
thousand  for  A,  and  one  in  three  hundred  ailiicn/cne  in  one 
million  for  fo  for  the  third  iteration.  The  program,  as 
presented,  steps  after  the  third  iteration  because  each  full 
iteration  requires  ^5  seconds  on  tae  HP-35  and  that  vas  the 
point  chosen  in  tne  precision/time  *r?.de-off  for  this 
particular  application.  This  fill  be  further  discussed  at 
the  end  of  this  chapter.  The  results  for  the  fifth 
iteration  (in  the  same  manner)  are  one  in  thirty 

theusand/one  in  three  thousand  for  Q,  one  in  one  hundred 
thousar.d/one  in  five  thousand  for  A,  and  one  in  five  hundred 
million/  one  in  twenty  million  for  fo. 

Precision  has  bean  discussed  since  it  is  independent  of 
the  equipment  utilized.  Accuracy,  on  the  ether  hand. 
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involves  net  only  the  precision  tc  which  any  given 
instrument  measures  is  wall  is  that  of  the  overall  system 
bur  also  the  accuracy  of  each  measurement  and  the  accuracy 
of  the  overall  system.  For  example,  this  system  will 
measure  the  given  value  of  the  temperature  to  within  more 
than  five  significant  figures  oetween  -90  degrees  and  150 
dearies  Celsius,  '/at,  it  is  only  accurate  to  within  six 
degrees  Celsius  (a  function  of  the  thermistor)  .  Thus,  there 
is  a  dramatic  difference  between  precision  and  accuracy. 

A.  CONCLUSIONS  AND  3 2  CONM  2NDATI 0N3 

The  ability  of  this  system  to  obtain  2,  A,  and  fo  far 
surpasses  that  of  any  manual  aetned.  In  addition,  it  allows 
the  operator  to  perform  additional  functions  as  the 
experiment  progresses,  does  not  make  transcript  tor./ 
interpretation  errors  as  is  ?o ssible/liheiy  with  manual 
recording  of  data,  gives  highly  reproducible  results 
independent  of  the  shill  of  the  operator,  and  allows  for  the 
resulting  data  to  be  to  be  easily  presented  in  a  variety  of 
formats  in  a  very  short  period  of  time. 

It  was  originally  intended  that  this  complete  system  be 
a  prototype  to  then  be  transferred  to  a  system  utilizing  the 
HP9836  rather  than  the  HP-35.  However,  the  equipment 
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arrived  too  late  to  be  incorporated.  One  major  advantage  of 
the  HP9836  is  that  it  has  an  eight-to-one  speed  advantage 
over  the  HP-85.  Speed  is  the  reason  that  only  three  rather 
than  five  iterations  of  'Ravine’  were  utilized.  Since  each 
mode  too’*  apprcxiaat  ely  135  seconds,  6.75  minutes  elapsed 
from  the  start  of  Measure',  lode  1,  tc  the  next  start  of 
'Measure',  Mode  1.  Obviously,  the  greater  the  number  of 
modes  the  greater  the  time  required.  This  time  diffference 
could  be  a  major  proolem  if  the  changing  external  parameter 
caused  the  resonance  of  the  moies  tc  shift  significantly  sc 
that  the  response  no  longer  fitted  within  the  chosen 
bandwidth.  One  method  of  overcoming  this  would  be  to 
incorporate  a  prediction  of  -he  next  location  based  upon  the 
change  of  the  parameter.  This  was  net  incorporated  in  this 
program  for  two  reasons.  One,  iz  was  never  necessary  in  the 
experiment  run;  and  twe,  the  shift  to  the  HP9336  with  its 
speed  advantage  woul 2  eliminate  this  lengthy  time  period  (it 
could  perform  twenty-four  modes  in  the  same  period  it 
presently  takes  for  three). 

Another  future  change  could  be  to  work  in  'Ravine'  for 
the  first  set  of  data  while  taking  the  measurements  for  the 
second  set  of  data.  This  would  be  a  better  utilization  of 
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computer  time  since  there  are  currently  periods  when  the 
computer  is  idle  (tfAIT  periods  and  time  constants). 

Another  improvement  would  be  to  replace  the  Icck-in 
analyzer  with  a  computer  controllable  model.  This  is,  of 
course,  a  trade-off  of  price  versus  cor. veaience/time  since 
there  is  approximately  a  two-to-one  ratio  in  equipment  cost. 

There  are  also  numerous  computer  programming  techniques 
which  could  shorten  portions  of  the  program  and/or 
associated  run  times.  One  obvious  change  would  be  to  make 
CRT  presentations  of  graphs  optional  or,  if  the  system  under 
test  was  already  well  known,  to  make  the  option  available  to 
enter  the  program  at  the  ‘NEW*  TAO'  point  thus  eliminating 
the  need  to  establish  the  desired  mode  data.  This  was  not 
done  in  the  presented  program  since  it  will  be  utilized  by  a 
variety  of  users  and  each  will  be  able/is  encouraged  to  make 
modifications  for  specific  applioa  tions. 
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A PPENDI X  A 


REFERENCE  GUIDE  TO  PROGRAM  LISTING 


Lins  Nos. 

1  -  251 

270  -  640 
650  -  1690 
1651  -  189  9 
1900  -  2676 
2679  -  4997 
4998  -  5320 


Ds  scription 

IN  PUT 

SEARCH 

SORT 

NE»  TAU 

MEASURE 

RAVINE 

TRACK 


53 


1  DIM  F h  {  20 G )  ,A  <300)  ,C(  ICO) 

2  CLEAR 

3  PPIN1ER  IS  2 

4  OUTPUT  722  ; "F4M6" 

5  CISP  "REMOVE  PROGRAM  TAPE  AND  INSTALL  CAiA.  1APE"  t  EEE 
? 

b  CISP  \vriEN  COMPLETED  ,  ENTEF  l" 

7  INPUT  Y 

6  If  **1  THEn  5  ELSE  CLEAF  t  SEEP  -  CC1C  2 

9  EKASE1APE 

10  CREATE  "DATA"  ,bbU,oo  s  REWIND 

11  CLEAR  v  LEEt  !  mtiai  sec  up 

13  SETTI ME  u  ,0 

21  DISP  "Enter  lcwer  tree" 

31  INPUT  fl  i  lcwer  tree 
41  CLEAR  «  EELP 
51  CISP  "Enter  upcer  rret," 
ol  Input  F2  !  upper  rrec. 

71  IP  F  2  <  =  F  1  THEN  CLEAF  is  SEEP  fe  CGTC  51 
61  1*1/  (  4*  (  (F2-F  1) /25b)  )  !  tin.e  constant 
Pi  CLEAR  is  SEEP 

101  DISP  "Set  tiir.e  ccnstant  on  5204  at  cr  siraller" 

111  DISP  "When  ccir.plete,  enter  value  set  (in  i'ilii-sec)" 
111  INPUT  T 1 

131  IF  T  1>T*  1GCC  THEN  CLEAR  5  SEEP  is  GOTO  10 1 

14  1  CUE AR  ts  EEEP 

1"1  CISF  "Enter  air.p  ct  cnv.jna  rrec  in  rrv  ( =<  3  50 C )  RMS" 
161  DISP  "CECIMAL  VALUES  ARE  NOT  PERMITTED" 

171  INPUT  A  !  airp  cr  rrec 

131  IF  A>  3  50  0  THEN  CLEAR  Is  BEEP  <s  GOTC  151 

16  2  FOR  N=1  TO  9 

18  3  A8(N)-A 

164  NEXT  U 

1S1  CLEAR  i a  EEEP 

201  CISP  "Enter  largest  value  ct  an.p  ever  cesirea  (.=<350 
Oir.v  RKS )  " 

211  INPUT  Al  !  n.ax  tuture  air.p 

221  IF  Al>  350u  THEN  CLEAR  ts  SEE P  (s  GOTO  201 

230  !  tnese  are  tne  values  tc  set  up  tne  3325A 

231  A$*vAL$  (A) 

232  A2$*"A«" 

23  3  A  3$= "ME" 

240  OUTPUT  717  ;A2$,A$,A3$ 

241  WAIT  200 

251  F3*CEIL  (  tF2-Fl)/256)  1  smallest  int  >=  bancwiotn 

252  !  SEARCH 

253  !  SEARCH 

254  !  SEARCH 


261  1  next  group  senes  free  tc  3  325A,  anc  gets  air.p 
341?  7A 
269  CLEAR 

276  DISP  "I  air  wcrxing  in  SEAKCH” 

271  FOR  A=G  TO  256 

2c I  F4 (K) =F1+K*?3 

2o2  1  c 4  ( K )  is  tree  sene 

290  F4$=0'AL$  (F4  (N)  ) 

291  H1S-"KZ" 

292  ii  2  $  *  "  F  R  " 

cUO  OuTPUT  /L7  ;H2$  »F4S  iiils> 

31u  rtAIT  4  *  T 1  + 1 0  C 
220  OUTPUT  709  ;"vT3" 

. -  cO  EATER  7u  9  ;  A  (A) 

2  5  u  A  E  XT  C. 

351  I  R5  AAL  t6  APE  SCfLii  FACTORS 

352  PEEP 

3oG  F  5*r  1- .  i  *  If  2*f  i ) 

37C  F  c=f 2+  .  1*  (  F  2.-F  i) 

: cC  GCLLAR  is  CLEAR 
-90  S  CA  ll  £5f£>0f~~«lfl«2 
4 0 u  XA/.Iu  o , 1CCC , F  1  ,  r  2 
4iu  TAXIS  f  l,.l  ,-.1  ,1.2 
411  FOR  A=G  TC  25o 
<*12  PER  OP 

413  PLOT  F  4  ( li )  ,  A  ( A ) 

4  In  ;,EXT  A 

s  15  ElSr  "  Ic  ICo  L£.;Ir L  rf  CCPx  Re.  TER  1  ,  EncE  E<;»riRR 
EcEi- 

416  IXPuT  x 

•ill  IF  1=2  THER  47  7 

4  io  GCLEAR  t  CLEAR  y  EE EP 

419  LISP  "EATER.  FIGURE  5  TO  EE  PRlATEu" 

<i20  IRPui  LE 

421  PLOTTER  IS  703 

422  PEA  1 

42  3  SCALE  t  l-.l*  (F2-F1)  ,F2  +  50G,-.. 2  ,1.  3 

424  XAAIS  u,l0Cu,Fl,F2 

425  TAXIS  Fl,.l,u,1.2 

426  PEA  2 

427  FoR  A*U  IC  250 
426  PEN uP 

429  PLOT  F<*(A)  ,A(A) 

4  30  A  EXT  A 

431  tERLP  y  PEA  1 

4  32  LEIR  C,SIN(9U)  is  PER  UP 

4  i  -  FOR  X=F  1  TC  F  2  STEP  2CuU  y  PEA UP 

4  34  MG  v  £  X  ,  - .  1  3 

4:5  LAEtL  vALS(X) 

4;o  i.uXT  X  A  PER  UP 


r  reir 
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4 4 C  LL IK  L 

441  FGF  i*U  PC  1.2  STEP  .1  t  PEIiUP 

442  CCvE  F 1-. U 5* (F2-F i)  ,¥ 

*■.**:  LAcEL  *A.LS(1) 

444  XEX'i  i  y  FEN  OP 

450  LLIP  C  y  PENUP 

4  51  i'.C  v  £  r  I-tjOGU  ,-.21 

45  2  LaeEL  " t FEwOENCY  IN  tio“ 

452  LL  IF  0  y  PENUP 
45  4  ;-’C  V E  FI,—  .  2 

455  LAE&L  “PLwI  Cr  bcLAl iv  S  A.*rLI  iwLc.  »;  (FULL 

SPECIE  Li-  )  " 

45b  LEIt  L,SIU(5C)  l-  PEi.Lt- 
457  KCvE  r  L- . Ob*  ( F2-r  1)  ,  .  : 

*oo  LAoEL  "PEL AT  I  vE  AnPLI'IUEu" 

-0  3  Ll IP  L  y  t-ENUP 
4  b  C  PC  VE  c  j  -f .  uCL  1 1  .  .. 

4b  1  LAE EL  "FICUFE  “;Ls 
*tL2  PENUP 
47b  PLOTTER  IS  1 
i  /  7  CEEAr.  -c  oLtt 

*t7o  LISP  "Lc  vcl  aesire  tw  rtri-r,  exj;  wit:,  cirierer.t  Cora 
retect" 

47  p  CISP  “U  sc  enter  1,  otherwise  enter  2" 

*!  b  U  I  Ii  P  L 1  I 

•» 5 0  GCwEAP  y  CLEAR  *  SEEP 
:i/u  Ir  i*i  liitN  li 

510  LISP  "Enter  cecisior.  ccir.t  ter  cCtiKtct  ( v  *C  1.2;" 

52o  INPUT  A 4  I  CECISICN  PCI M 
5:0  CLEAf-  v  LtEP 

54o  LISP  “Enter  cull-scale  cer.ai  civ  it/  setting  iron:  cne 
52l4  ir.  vcits" 

550  INPUT  A 5  !  SENS II I v 1 1*  SET! INC 

5c  0  PRINT  "An.plituce  in  vclta"  ;  “  frecuency" 

570  PRINT 

5oO  FCR  N*0  TC  25o 
5S0  IP  A  ( N ) < A 4  THEN  o:0 
bOL  A  (N )  =A5  *A  ( ti ) 

olO  r PI M  USING  bit  ;  A(I\)  ,t  4(N) 

o2u  I FACE  IX  , L . LLLCCLLL , 1  OX .CDCLLC .CD 

bit  NEXT  N 

631  PRIM  USING  6  22 

o:2  I  FAGS  1/ 

t40  CLEAR  y  tEEP 

641  1  SCAT 

642  !  5CFT 

643  !  SOFT 

650  LISP  "hew  irany  ircces  c*l  you  cesire  tc  tracx  (t-Ax  CF 

9)  “ 
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6  6  6  INPUT  t.  !  NUMBER  CF  NODES 
b  7  0  tCF  A *  1  TC  K 

bbC  CLEAR  t  eEEP 

650  DISP  " rvnat  is  center  tree  tor  ircce  ";f>. 

7 00  INPUT  t1.  ( K )  !  CENTER  FREW  " 

7 lU  CLEAR  «  PEEP 

7  20  DISP  " An a c  is  rrecwietn  ter  n.cce  ";i< 

7  30  DISP  “ RUST  EE  GREATER  THAN  "  ;2*F3  !  F.2  IS  £.. 

746  INPUT  C(A)  !  FREu~ WIDTH 
7  5  0  RE  AT  A 
7ot  CLEAR 

761  !  cren.  77u  tc  560  is  tne  tirst  rcutji.  try  rcr  ir.easuri 

ne  (SC  FI) 

776  FC  T  L=i  TC  .V 

77L  lISP  "I  an.  wcrxing  in  SCRT  ter  nece  "  ;L 

7c 0  FOP  A=0  TC  lob 

7  Sc  F  4  (K  )  ■».;  ( L)  -C  {  L )  /  2+i  *C  (L)  /  let 

oCu  t 46*vAL5 ( F 4 ( 1 ) ) 

bio  CUT  PUT  7  L  7  ; H 2S  ,  t  b  ,ri  1 S 

6  26  rtAIT  ik^Tl+loo 

o  20  OUTPUT  765  vi  3" 

o  2  1  rtAIT  5u 

b40  ENTER  70S  ;  A(N) 

E5o  NEAT  u 

651  LEEP 

co  1/  GCLEAK  i“  CLEAR 

b76  SCALE  t  4  ( u)  -  .  1*0  <L)  ,F 4  (  it  t)  +  .  1*G  { L )  ,- .  i  ,  i .  2 

obO  AAAIS  v,,C(L)/ic,r<*(0),cMiOO) 

bSO  YAAIS  F 1 < u)  ,  .1  ,6  ,1. 2 

i'  o  6  rGR  A =  6  a  U  lot 

Sit  PEaUP 

S2C  PLC'l  F  4  ( A )  , A  { N ) 

5  30  A EAT  A 

c*40  >.CVE  F  4  (  0 )  1 

54  1  LAcEL  'VC CE  ";L;"  "  ;  F  <.  ( C )  ;  "  TC  '';F4(  106) 

“42  LlSt  IF  YOU  rtAAT  rt  CCPi  EA  1EF»  1  f  Et>SE  tATEF.  2 11  f  >  E£ 
tP 

54  3  INPUT  Y 

54 4  IF  T=2  THEN  So7 

54  5  GCLEAR  is  CLEAR  is  oLEP 

S4o  DISP  “Eater  FIGURE  *  TC  EE  PRINTED" 

54  7  INPUT  LS 

55  4  P LC ITER  IS  7U5 

5 55  PEA  1 

S5o  SCALE  F4  <0).-.2*G(L)  ,F  4  ( It C)  +  .  1*0  ( L)  ,-.3,1.3 
»57  /.AaIS  6  ,C(L) /lu  ,F4  (u)  ,F4  (  106) 

556  *AAIS  r 4  (u)  ,  .1  ,6  ,1. 2 
555  PEN  2 

5tt  FOR  A  =0  iC  10 o 
»oi  PEA6P 
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Sfcl  PLOT  f  4(N)  ,A (h) 

90  3  NEXT  K 

964  FE*.. LF  (y  FEN  1 

9o  5  LC  IF  t,SIN(90) 

966  F OF  X=F<,(0)  TC  F4(10t)  STEF  C(L)/1C  t  FENLi' 

9  o  7  AC  v  E  X  ,  —  .  1  3 
ioo  Lr  tEL  VALv(X) 

969  NEXT  X 

97t  LEIF  t  t  FENLP 

9*1  FCF  Y*u  TC  1.2  STiif  .  i  i  FcNLF 
9  / 1  r.C  vt  c  4  (  G)  •,  ct  *C  ( L )  i  i 
9-7  3  LABEL  vALE ( Y) 

9/4  next  i 

91 7  5  LCir  t  i;  FENCE 

9/0  AC  vE  t4(G)4.l*o  (L)  ,-.11 

9  77  LAbEL  " F FEuoEXCi  IX  rii." 

9/1  L  n  I  f  L  FE'iLF 
9  79  ,VC  V E  c  4  (  U )  -  .  1  *C  ( L )  ,  -  .  3 

9oC  LAbEL  "  F  LCT  Cc  i-  lLAI  I  v  2  A.  FLIT  LEE  vS  FFEwL  Li’Ci  ,  .-/CEE 

»  .  r 

/  -oi 

9  cl  L£IF  L,  SIX (90)  ^  FEi.LF 
96  1  ^  u  vE  f  4  (  o ).-.  1 7*C  (L )  ,  .  3 
9- c  3  LAtEL  "ILLATIVE  ALrLIILLE  " 

9  6  4  LEIF  L  t  PEKLP 

9c3  i-.CvE  t'4  (  0)  +  .1*0  (L)  ,1.2 

9  co  LAE  EL  "fc'IGLPE  "  jLS  9  rt.'.Cr  o-  P  LETT  EF  IE  1 
9c7  CLEAR  £  oL.Lt 

900  u  X  S  t  OC  jf  a  c.li  t  Co  ^  Si.yonl  J"iO»  .  I  £  wC  orccc  X, 

ci£t  i 

3o9  Ii\Pu.'I  ic  CLoArv  ELEF 
990  Ir  9  =  2  'In EX  120o 

iooo  EISF  "Cnsr.^e  seits  i  c  l  v  i  cy  ?  Enter  i  £  cc  yes,  2  ror  r.c 

II 

101c  lx, ALT  it  CLt,AK  s  c££F 
io20  It  i  =  2  Xii L X  lobe 
lu  :o  CLEAF.  t  3LEF 

lo40  uISF  "Enter  nev,  sensitivity  setting" 

IE  50  I.,  Pul  At 
loot  CLEAR  t  EEEF 

lo7C  D ISP  "Cnar^e  £recv*ietn?  Enter  L  rcr  yes,  2  rcr  no" 

loot  IXroT  it  CLEAF  ^  EEEF 

io9E  IF  i  =  2  Ir. EX  1120 

lloo  C IS r  "Enter  now  rrecwiccn" 

lilt  I L« PuT  C(L) 

112o  CLEAR  e  EEEF 

ii30  C  ISP  "Cnantje  center  tree?  Enter  1  rcr  yes,  2  rcr  r.c 
*• 

1140  1 1. ALT  Yt  CLEAF  i  BEEF 

1130  If  i =2  THEN  llct 

1 1  b  t  E  IS  t-  "Enter  new  center  tree" 
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1170  II'.r'LI  i-  (  L ) 
llot  GCLtAr.  t  C-LcAr- 
1  Ho  GC1C  771 

11*1  l  trcit  L2Uo  tc  Lcco  is  tr.e  tiret  tr>  at 
ter  tree  arc.  tr.e  s* 

120o  Clitftf  t  CCL£rF. 

12o  1  CIS*  "I  tji  i»cr\ii‘.'j  ct.  tree,  ar.c  v  j-^-r  £C 

121u  ri  ,cl  ,E4  ,ti  1  ,r.2  ,r:  7  ,ti  t*o 

1  <.  2  o  22  ,£2*50 

1  4.  C  0  f  «  *  H  *  1  A  C  i  C 

Ucc  If  A(.N)<A(ij)  i:2i-  Lee 

122  1  If  A  (.1)  »A  (k)  iilti.  in  2  cLSC  12  5  u 

lit  2  rrl.i  l  A  (L) 

1 2  5  o  a  o  ( i. )  *A  ( i\ ) 

12uO  r  o  <L)  *Ft  (!'  J 

1  *.  y  o  ii*i» 

12oc  1.2X1  1* 

lit 0  A7(L)*Ao(L)/Sv.f(2) 

1:00  to?  o.  =  C  IC  i. 

121C  If  A  7  (  L)  *A  (2  )  inil  Hit 
1  :*.o  If  A  /  ( L )  < A  ( I. )  1 1: c !•  1  :  /  c 

13  3c  Ir  A  ( N )  <  t  j  It? 21.  i<IDi 

1  : 2 0  ci*A(rt) 

1.350  ijl*L 

13  co  ssCfC  12  50 

1  3  7  0  It  A  ( L )  >  1 2  fr.ee* 

13cc  2 2* A  (;••) 

13  50  « i  2*ti 
1400  Gv.  iC  It  CO 
It  10  C  L  ,  E2*A  { 1. ) 
it^o  HI  ,f.2~ W 
1230  F7*c4(m2) 

14 tO  GClC  147o 

It  50  wEXT  i-i 

1451  a»{A7(L)-£1)/ (fii-cl) 

It 6  o  F  7*X*  (f  4»  (is 2 )  -c  4  (II 1 )  )  +  c'4  (.i  L) 
l47o  fGF  N*h  IC  100 
14 oo  Ir  A7(L)*A(1.)  ’I tit L  15  CO 
1 1 5- o  IF  A  7  ( L)  >A  ( is )  it iE'h  1 5  1 0 
1500  IF  A(N)>c3  FUEL  15  ;o 
15  lo  c3*A(l») 

152u  ti3*N 
1530  GCTC  1020 
1 5 4C  IF  A (N ) <*B4  IKEw  157C 
15  50  £4*A(N) 

15C0  H4*to 
1570  GCTC  lo2C 
15cU  £  3  ,t4*A (N ) 

1550  ti3  ,H4*N 
loCO  Fo*F4(ril) 


161c  GCTC  16  40 

16*0  A  L.  A  1  i\ 

lb 21  X® ( A  7 ( L ) —  E 4 ) / (fi  :  — E 4 ) 

16  3  0  F&=- (X* (F4 (H4  ).~F  h  (til)  )  )  +F  4  ( h  4 ) 

16'4C  u  (L)  =  r  fc  (L)  /  ( F  U.-C  7  ) 
lb  So  PR  IN T  'VOCE  ";L; 

iboo  PRIM  "Ct’ViEF  FP£<*  IS  “;F6(L);"  AND  A:P  I£  ";A6<L)* 
A  5 

It  7  0  PRINT  '\  IS  "  ;  U  ( L ) 

16  71  CLEAR 

i  b  7  *  F  R I  n  T  uS  I  n  C  lb  >  3 

16  7:  IMAGE  5/ 
loco  EXT  C 

1  o9  0  CLEAR  GCLEAP 

1651  CISP  "I  an  wcr;;ir.g  cr.  tir.e  ccr.scc.nc" 

lb5S  !  Frcrr  17o0  tc  17^0  ere  largest  Tire  Ccr.scai.c  is  sc 

ur.c 

17  00  FUR  L—  1  'IG  M 

1701  T (Lj  =G (L)/ ( P I  *r  o ( L )  } 

1702  NEXT  L 
1  705  T-l  (  i ) 

1710  FCk  L=2  TC  .«• 

1711  IF  T>1 (L)  'I PEN  1720 

1712  T=T (L; 

1  /  *  0  i.EXT  1 

L 7 50  c.EbP  t  CLEAR 

17  55  PFIN'i  “TIME  CONST  >=  "  ;'i*Ku 

17oo  CISP  "Sec  ana  enter  r.ew  tire  constant  (in  rillc-sec 
J" 

1770  CISP  "FUST  EE  GAtATEF  THAN  ";1;H  sec" 

1760  INPUT  1 
17  5  0  T  1  =  E 

1751  CLEAR 

1752  PFINT  "MEN  THE  CONST  IS  ";'Il 

1753  PRIM  USING  1754 
17  54  I  PAGE  4/ 

1755  1  MEASURE 
175o  1  MEASURE 
1757  1  MEASURE 

1756  1  MEASURE 

1755  CISP  "I  air.  wcrxir.g  cn  scar.carc  cev  ,  #  cr  eta,  &  var 
iar.ee  " 

1800  !  Frcir.  ltOO  tc  1^24  is  tr.e  calc  cr  tr.e  noise,  st.  c 
ev.(  ar.c  tr.e  r.oirxer  cr  pts.  usea. 

IsOl  G,a=0 

loO  2  F*  ( F 6  (  2 )  -F 6  (  1 )  )  /  4  +F  6  (  1) 

160:  F4$  =  VAL5(F.) 

Ifc04  OUTPUT  717  ;h2S  ,Fh$  ,iil5 
lc.0  5  ..AIT  Tl*  12+ioG 
160o  OUTPUT  705  ;"vTi" 
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ltd 7  WAIT  lOO'U 
1 6  0  6  FOR  a  =  1  I C  1  o  U 
loOS  CM T ER  70S  ;  A ( K ) 

Itlo  v\*W+A(N) 
ltll  NEXT  N 

1612  *»=*»/  luC 

1613  FGF  N  =  1  10  10 0 
1  c  1 4  S*(A(w)-vv)*2 
1615  G=G+S 

181b  xEAT  N 

lo!7  G*A5*Ss<E(o/S:>) 

lb  16  iism*A; 

lc  IS  PRIM  "Tne  rear  ie.  " ;  ,v 
lclu  PRINT  "Trie  star.caru  cev  ic  ’*  ;G 
lcll  PP 12.  i  uSIaG  1621 
loll  Ii-AGii  i/ 

1615  FOR  L=  1  1C  i' 

16  30  I (L) »1.&2*SwR (A 6 (L) *A5/G) 

1635  PRIM  "PC.  OF  POINTS  rCr  ..CEE  ";L;"  I J  ";IiLJ 

1840  IF  I(L)<lou  TiitN  lc30 

1  c  4  5  I  ( L )  *  1  o  0 

165o  IF  1(C) >50  THEN  looO 

1655  I ( L) -50 

1C6C  I (L) = I P ( I  (L)  ) 

lc6l  I  ( L) *24 

lc  b  4  PRINT  "NO.  CiiCSEN  18  "  ;I(L) 

lc65  PRINT 

lc co  NEXT  L 

lot/  CLEAR  4r  PEEP 

i86  0  PRINT  OSImG  lObl 

1  o  o  1  I  FACE  4/ 

locS  PRICIER  IS  7 U  1  ,132 

loo 7  ClM  F $.  (  2 0  ] 

166c  FS  =  "c-G  FA v I  ME " 

lcSC  PRIM  T  USING  lbSl  ;  "III'.E"  ,  "TEf  P "  ,  "  Pf  ES "  ,  ".V 
REw"  ,  "APR  ,'V  ,  "SNR"  ."C-AMP"  ,'*5  FT"  ,Fo 

1691  IMAGE  4A  ,  5X  ,  4A  ,6  a  ,5A.  ,  4X  ,  3A  ,6X  ,  7A  ,  5X  ,4A  ,  L-4X  ,  1A  ,5 
,5X  ,  5A  ,4X  ,4A  ,6X  ,  luA 

1692  FEINT  USING  IeS3  ;  "(SEC),,,"i  C  HZ )*’,"( V  rn 

"  ( it  V )  " 

lc93  IMAGE  5A  ,4X  ,5A  ,25X  ,  4A  ,  1CX  ,  cA  ,  30X  ,4A  ,/ 

169  6  PRIMER  IS  2 
1699  CLEAR 

l9oo  !  Frcir.  I960  tc  4S6O  is  t r.e  calculation  icr  tne 

ore  inducing  center  tree  ar.c.  w 

1  SC  1  15*1 

190  2  FCR  L*  1  IC  K 

1903  CISP  "I  M  IN  MEASURE  FOR  MCCE  ";L 
1900  r  1  =  F6  (L) -Fb  (L) /w  (C) 

1SC7  t  2*Fb  (L)+Fb  (L) /<»  (L) 


C.  F 
X  ,  3A 
£)  "  . 

•Teas 
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1506  u=  (F 2-F 1)/ I  (L) 

1505  v  =  {  F  2-F  1)  /2 
15  10  ri=U 

IS  1 1  If  Ao(L)>.j  Ih£iv  1  if  1 7 
1912  A <5  ( L)  ®2*Ab  (L) 

1513  Ao {L)=2*A6 <L) 

1914  IF  Ao  (  L)  <A  1  THE:.  1911 

1915  A8 (L) *A 1 


1  5  16 

GCTC  1525 

1517 

IF  Ao ( L) <  .  55  TtiEF 

1525 

IFlc 

A 6 (L) =Ac ( L) / 2 

iFiy 

Ao(L)=Ao(L)/2 

1  f  20 

GCTC  1517 

1525 

A$  —  vAL5(Ac(L)  ) 

1526 

GolFoT  717  ;A25,AF 

,A 5 

1527 

v»A  IT  4l0 

1  f2« 

Ei.  'i  EF  7  22  ;  T  6 

1525 

FCF  i.  —  G  iC  I  ( L ) 

1 5  :u 

F  4 (N ) =r  o ( L) +l*C- v 

15:1 

F 4$-VA  LS  (r  4  (F  )  ) 

15  s* 

h 1$=" h 2" 

1  V  5  3 

h'  2  $  *  "  F  t\  *' 

19*54 

CLTFL'I  717  ;  f:2  5  ,F  4 

5  ,uU 

1  r  :  5 

aA  I T  1 2  *  T  1  f  1 0  o 

15  lb 

CoTFOT  705  ?  "  v‘i  3 " 

1  5  17 

EivTE B  705  ;  ( :•  } 

l  5  :  b 

-j  L  a  I  N 

1515 

c.v'I  n.r<  /  *.  4.  i  i 

If  4  0 

FC r  ..  =  1  iC  I  (L) 

154  1 

IF  A  (X  )  <A  (tt)  T;iEx 

1  i  6  0 

154  2 

IF  rt(I.’)*A(li)  Into 

1 5  4  1 

154  1 

FFINT  A(N) 

1550 

H«N 

156  o 

NEXT  X 

1570 

X*A  Itt ) 

19 tC  T*(A  (K  +  l)  -A(H-l)  )/2 
1950  2*  (A  (ii+i)  +A  (H-  1) -2*X)/ 2 
2uoU  c 6 ( L) »- ( i/ ( 2* o ) ) 

2olo  A6 (L) *X+1 *F6 (L) +o*F  t  (L) *2 
2  0  2  0  r  o  ( L  )  =F  4  ( A )  +r  5  ( L )  *  ( F  **  ( )  -f  4  ( H  - 1 )  ) 

2o  40  ,c  l*5o 

2290  A 7  ( L )  «At>(L)  /S*F(2) 

21oo  FCF  ::*o  iC  H 

2  lie  IF  A7(L)*A(M  THE:.  2<*lC 
2120  IF  A7(L)<A(N)  lttEi^  2;7o 
21:0  IF  A(iyJ<tl  T  tit  Li  2*t5G 
2  1  <*  o  ol=A(N) 

2150  ril*N 
2It.o  GCTC  24  50 
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2  37G  IF  rt(N)>02  liiEN  2*00 

1  -  b  C  t  i. — A  ( i  v ) 

2  3St  H2=N 
iiCu  GCI C  2450 
24  It  B1  ,n2=A  (N) 

24  <.0  Hi  ,ii2=N 

24  30  F7  =  F4 (n2  ) 

*:44G  vjGIC  <.470 

24  50  NEXT  N 

ii  S  1  a=(A  7  ( L)  -  l  j)  /  (tl'c  ;  ) 

*.iuL  r  7=  X  *  ( F  4  ( i;  2 )  -F  4  ( i:  1 )  )  +  c  L ) 

<•4/0  irCF.  i.  =H  iC  I(L) 
iioi/  I c  A.  i  (L)  —A.  ( i\ )  Tb.EN  l5cI 
Itit  IF  AT  ( L)  >rt  {  >. )  'ir,u.  25*,b 
<5tt  IF  A(i'j)>c;  'Ir.Li*  25:0 

25  It  1.3=;,  (N) 

2  5  2 1  h  3  =  f . 

«. 5  : 0  CCiC  i. b i o 

i. 5  ^ 0  IF  A  ( .» )  <=  0  -t  T hEi.  «. i  7 1 

25  50  E<i=A.  (:») 

1  ;  o  t  r:  4  =  2 

2570)  GOTO  2olC 
2500  c  3  ,c4=A.  (N ) 

25  90  h3  ,h4  =  N 

<.  o  b  t  F  c  *  F  4  (  n  3 ) 

2o  It  GU'IC  2o*ti/ 

4.  b  4  O  L  X  T 

2021  a  =  ( A.  7  IL)”ti)/(t:“c<i) 

<-6  30  t  os*(n*  ( F  4  ( H*»  ).~r  ft  (n  3)  )  )  +F  4 ( n  < ) 

2b4t  UL)*rO  |L)/(ta-r  7) 

2b  7  :  PRIMER  IS  7t  ]  ,132 
2fc7b  CLEAF 

26  7  7  ! 

26  70  ! 

267b  !  vAFY  * 

2  bob  FOF  0=1  TC  3 
2  b  a  5  FCF  K=1  TC  3 

+  otb  GI3P  "I  AM  IN  FA  VINE  FCF  v  OF  wCEE  " ;  L 
269 1  o(l)='*(L) 

27 tu  0 ( 2) =1 . GC  5*o (G) 

2710  0  (  3)  =.  9b5*UG) 

2711  F*A6{L) 

2712  0=( 1/0 (X) ) *2 
2720  E ( K) =t 

2730  FCR  N* 1  TO  I (L) 

274  0  C (N) -R/ (0  (K)  *S*F  (  (F4 ( M) /F 6 ( L) -F 6 (L ) /F4 ( N ) ) *2  +  0 )  ) 

2750  D  =  (A  (N)  -C  (N)  )  "2 

2760  E(K)=E  (K)+D 

2770  NEXT  N 

2775  CLEAF 
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2760  NEXT  K 

2813  v  =  -(  {£  (2)  -E  (  3)  )/(2*  (E  (  2)  +E  (  3)  -2*£  ( 1)  )  )  ) 

2614  Q{L)=J  (l)+v*.C05*J  (1) 

287b  ! 

2877  I 
2676  ! 

2679  !  VARY  AM> 

2660  CLEAR 

2661  FCR  K« 1  TO  3 

2669  LISP  "I  Ax-1  IN  RAvIRE  FOR  Al-PLITUDE  OF  RCC  £  ";L 
2690  0 { 1 ) *A6 (L) 

2900  J <2)=1.0t2*A6 (L) 

2910  0 ( 3) =.99c*Ao (L) 

2911  Li*  ( l/w  (L)  )  -2 
2920  £ ( K) =0 

2  9  30  FCR  l\=1  TC  I  ( L ) 

2940  C  (N  )  *J  (K )  /  (w  (L)  *S  wR  (  (f4  ( I. )  /  t  o  ( L )  -Fc  (L)/F4  (1. )  )  2+U)) 

2950  D=  (A  (K)  -C  (N )  )  '2 
2960  £ (K) =E  (K) +C 

29  70  NEXT  xN 
2975  CLEAR 
2980  NEXT  K 

3020  V  =  ~(  (E  (2)  -E  (3)  )/(2*(E(2)+E{3)-2*E(l)))) 

30  30  A6(L)*J  «1)+V*.0u2*u(l) 

317  6  ! 

3177  1 
317  6  ! 

3179  !  VARY  FREQ 

3160  CLEAR 

316  1  FCR  rv=  1  TC  3 

:  lo  5  CIS?  "I  Ax-1  IN  RAvINE  FCR  FREQUENCY  OF  NCCE  "  ;L 
3190  u ( 1 ) *F6 (L) 

3200  j(2)*.005*F6(L)/s«(L)+F6(L) 

32  lo  J ( 3) *-(. U0  5*Fb (L)/Q(L) )+F6(L) 

3211  f.*Ao  (L) 

3212  o=  (1/Q (L) ) "2 
3220  E ( R) *0 

:2  30  FCR  x\  =  l  TC  I  ( L) 

3240  C  (N)  *R/  (C  «L)  *SQR  (  (F  4  ( N) /J  ( M.-J  (  R ) /F4  (N  )  )  '>2+o  )  ) 

3250  D=  < A ( N) -C (N )  )  '2 
3260  E ( K) *E  (K)+C 
3270  NEXT  N 
3275  CLEAR 
32SO  NEXT  K 

3  30  5  v=-(  (E  (2).-E  (3)  ) /  (  2*  (£  (  2)  +£  (  3)  - 2*E  ( lj  )  )  ) 

3306  c  t ( L ) *0 |l)+v*. 005*0 ( 1) /u  (L ) 

3:07  NEXT  t 
3  306  ! 

3:09  ! 
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3  310  I 

3  311  1  LAST  SiiCT  RAv'IXE  t, 

3  312  CLEAR 

3314  LIST  "I  AAi  lii  LAST  SiiCi  K" 

3  3 lo  u  =Ab (L ) 

3  3  17  L=  (  i/'v  (  L)  )  a2 
3  316  c  =  u 

3  316  FOR  W»1  1'C  I  (L) 

3  326  C (N) =J/ (w (L) *SwE ( (r  4  ( X )  /F  6  (  L )  -r  6  (L  )  /  E  4  (  R  ) )  “2  +  0 )  ) 

3  32  1  L=(A('.)-C  (X  i  )  *2 
;.2i.  E-S  +  L 

;  :Zo  u  —u  *A  5 
;  -  <-  /  T  =  o  /  vS 
3  32c  L  12'  06  [  1  32  ] 
j  .3  5  c  a w^IGi.7  1  TO  11  LcilA'* 

3  3  o  0  14=11  ;v.  2 
3  37  1  T  7  =  (  T  5+ T  c  )  /  2 

«  i>  t"w  =  0 

3  3  5  0  ThliTl?  1,15  ;  T4,l7(A3/L,F6(L),.J,o(L),l,Ac(L),I5,L 
342u  ASSICR*  1  iC  * 

—  4  0  1  Oi  *3  “  -b  |  ^  /.  f  •  U  i  T  A  A  i  -wL  |  |  ihL  (  w  f  TC  ^  t  i  vy  ^  4  /  /  O  •  4  L  C. 

I  0  A  /  3  L  •  40  i  3<\  i  3  LO  1-^  I  4L  i  i  /.  i  -L  |  C  A  fti  *:v  t 

3 40 2  RFIM’l  oE  IRC  os  ;  1  4  ,1  7  ,R  5  ,L  ,r  t  i L)  ,o  ,  v  ( L)  , x  ,A 6  (  L )  ,  1  3 

.ci 

l4io  15=15  +  1 

3“2o  r  •,  s>  =  v  A.  L  6  (roll) 

3  •=  3  o  A  6  =  v  A  L  v  i  A  c  (  l.  )  ; 

3  4  40  Col  ROT  7  j.  7  ;Aos,A6/A.  35 
3  4  ;  o  /.Ax  1  2  o C 

34oo  o/olRoT  717  ;i-.2s  ,s  4;  ,nis 

3  +  7o  a  A II  li*lo+loo 

34to  ColRoT  7o5  ; " v T  3 ” 

34  60  ERi'lER  705  ;  Cb(L) 

4  3j76  rF.IX'IEK  IS  o 
4 6 co  CLEAR 

4  6  c 1  NEXT  L 
4  5 1 2  R  R I N  T  E  R  IS  /  0  1 
4 So  3  RR.IiVi  LSIhOj  4  5b 4 
4  S  6  4  I  i-.A  0  E  / 

4Sb3  REIXTEE  IS  2 

4  56  5  ! 

455b  I 
4657  I 

45  6  6  I  TEA.Ck 

5ooO  I6=o 

3olC  FOR  L=  1  TC  >v 

Soil  CISR  "I  IFACRIXO  i-GLE  '' ;  L 
5o2U  f  4  i>  =v  A  LS  (  r  b  I  L  )  ) 

So  3o  A.5=vAL5  (At  (  L)  ) 


80 


in  in  c/'  tTi  u*  4 r  4fi  O'  tn  ( n  <n  u>  in  cr*  <»i  i.n  in  in  m  ui  m 


uulFi'i'  / 1  /  ;Ai*|A^  ((i 

a A 1 1  2 Ct 

Gt'iPta  717  ;  h2  v  ,  f  4  Z  ,r.  1  2 
hA.II  i  1*  16t  lo  u 
GclAul  7 ;  “  v  i  3" 

Li/lch  7tS  ;  AS 

IF  A 9  <.=  .  2*C o  ( L )  2i iti-j  tilt 

IF  A2>=1.  i*Cfc  (L)  X HEI.  tilt' 

GCTC  52*»5 

X!»«l 

IF  a  2  >  .  3  i  r.  £  2 .  5  2 1 1 

Md  (L)  =  2*A  t  (  L  ) 

A  t  =  2  *A  t 

lx*  A  c  (  l*  )  *C  A  x  i !_ A  2  1  x.  c 
A  o  (  L )  =A  1 

IF  A  b  <-.  ii  IriEG  5x^5 
A  d  (  L )  —A  o  (  xj  )  /  x 
A  t=A  >/ 2 
GOT C  5  x 1 d 


'J  I  L»  w 

5  :t  1 


Mt/ll  L 

ui-. 'I El  7  22  ;  'io 
IF  A c 5 ( 1 b -  1  7 )  >  = . 
IF  1  F <  i  l i.G.v  5u  . 
AFI.ViLF  IS  x 
14-1 1  .■  E 


2 x  u  1  o  njrLL. 


■  E  n  x  x  x  . 


oiii 
3  t  d  0 


G  u  xX  1  3  t  x 
GLSAF  -i  GCLcAK 
EISA  "  'i  h  E  tl.L  " 


APPENDTX  B 


A  sample  program  is  lists!  for  obtaining  information 
from  the  data  taps  aid  utilizing  it  in  a  figure  or  plot. 
The  key  spot  in  this  program  is  at  line  31,  the  reading  of 
data  from  the  tape.  There  are  a  maximum  of  850  records  cn 
the  tape.  In  this  case,  15  is  indicating  which  record  is  to 
be  accessed.  Each  of  these  records  car.  be  accessed  in  a 
random  manner.  Only  the  actual  value  is  stored,  not  the 
specific  variable  than  represents  it.  Thus  for  example,  if 
in  the  original  storage  program  the  values  F  =  1,  T  *  3, 


G  =  4, 

and  B 

=  7 

were 

stored. 

only 

th 

e  actual  values 

w  ouli 

be  st or 

ed  in 

the 

or  der 

given 

{1,3,4 

.1 

).  tfher.  reading 

t  hese 

ValU9  3 

from  t 

'U 

(1) 

r n  e  1st 

value 

m 

ight  be  called  A, 

the 

second  Z,  etc.  The  important  point  is  that  if  the  4th  entry 
in  the  record  is  desired,  four  values  must  be  read  into 
memory  (read  n  values  to  get  ath  value)  and  whatever  the 
reading  program  calls  the  entry,  is  what  it  now  must  be 
referred  to  within  the  program. 


1  C LEAK  f  DEEP  ii  GCLEA.P 

"uCm  MANY  t-0 1  NTS  ARE  BEING  Ev  ALLA  ill,?  tfi'i  EE  NU.NEE 

•h 

3  INPU1  E 

4  CLEAR 

5  CISP  "  LN TER  in£  LOWEST  X Ei-iP  CESIKEL  CN  PLCT " 
o  INPUT  El 

7  CLEA.R  *  EEEP 

o  LISP  “ENTER  1HE  hIGHEST  TEMP  DESIRED  CN  PLOT" 

9  INPUT  £2 
lb  CLEAR  A  EEEP 

17  CISP  "ENTEF  THE  NODE  f  DESIRED" 

16  INPUT  G 
19  CLEAR 

21  LISP  "ENTER  THE  FIGURE  ?  LEE  I  PEL  TC  EE  PRINTED" 

22  INPUT  E9 

2  -  CLEAx* 

27  DIM  1 1  (bSl)  ,12 (B51) 

J.B  GCLEAR  cJ  CLEAR 
2  9  ASSIGN?  i  1C  "DATA" 

3C  FOR  15=1  TC  b 

31  READ?  1,15  ;  1  4  ,T7  ,P  5  ,L  ,F  6 

32  IF  L=G  THEN  4U  ELSE  50 
4 U  Tl(  I  5 )  =  (F6/G)  *  2/ 1G  ~b 

4 1  T  2 ( I  5 ) =17  +  273. lb 

50  NEXT  15 

51  FLCTTEF  IS  705 

52  PEN  1 

oC  GCLEAR  iy  CLEAR 
61  £7 “FLOOR  (b  1+27  3.  lb) 
o2  co=CEIL (£2+27 3 . 16 ) 

7G  SCALE  E 7- lo  , c6+ 5  ,  3 5  , 5 1 
cU  XAXIS  40 , 1G  ,E7  ,££ 

50  YAXIS  67  ,2  ,4G  ,50 

51  LEIF  0, SIN (9C)  fe  PEN UP 

92  FCR  X=6 7  TC  E8  STEP  10  S  PENUP 

53  MOVE  X  ,37. 5 
94  LABEL  VAL$(X) 

S  5  NEXT  X 

96  LDIR  0 

97  FCR  Y“4G  TC  50  STEP  2  £  PENUP 
96  MCVE  E7-4  ,Y 

99  LABEL  VAL$(Y) 

100  NEXT  Y 

101  LDIR  0  £  PENUP 

102  MOVE  67+1G ,36 

103  LABEL  "TEMPERATURE  (K)" 

104  LD IF  0  , S IN ( 90 )  £  PENUP 

105  MCVE  £7-6,37 

106  LABEL  " (FREW/M#) “2  IN  KHD^l" 
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jJLt 


I  U  ltd  A I  <1111  i  l  in 


] 


107  LDIR  C  4*  PEN  UP 
1  u  8  iMC  v  E  £7-8,35 

109  LABEL  " PLC1  CF  FEEL; “2  VS  TEKP  ,i-!CDE  "  ;G 

110  LCIF<  C  t  FEN  UP 

111  iv.OV£  87+5,50 

112  LABEL  "FIGURE  ";85 

149  FEN  2 

150  FCF  f,  =  0  10  E/2-1 

151  P.=G+  3*N 

152  PEN OF 

152  PLOT  12  (F)  ,T1  (P) 
lS't  NEXT  N 
155  PLNoP 

lcG  CLEAR  CCLEAR 
lfcl  RE wIND 

io2  DISP  “rtLL  l»-/N£  ,  1IR.E  IC  CC  CN  " 

1/0  ti*  N  D 
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